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CN ' Abstract 

H ■ 

Links between supersymmetric classical and quantum mechanics are explored. Diagrammatic 
. representations for ^-expansions of norms of ground states are provided. The WKB spectra of 

, supersymmetric non harmonic oscillators are found. 

<n ; 

^ ; 1 Introduction 

m ; 

In this essay, written to commemorate the sixtieth birthday of J. Carihena, we discuss several clemen- 
CO ■ tary issues in one-dimensional supersymmetric quantum mechanics. The role of the Riccati equation 

O ' 

\Q ' in this framework has been thoroughly analyzed by Carihena and collaborators at the highest level 
| of mathematical rigor by approaching topics such as the factorization method or shape invariance 
■i— > | from a group-theoretical point of view, see P, [2] and [3]. Our purpose here is to approach these 
Q-i, matters from a rather physical point of view. To construct a supersymmetric quantum mechanical 
system starting from a physical potential energy we shall be led to deal with the Hamilton-Jacobi or 
■ the Poisson equations, although in both cases there is an associated Riccati equation. We shall focus 
^ ■ on studying the relationship between supersymmetric classical and quantum mechanical systems, fol- 
lowing the standard References [3] and [3] and the more recent Lectures of A. Wipf [H]. In particular, 
models where supersymmetry is unbroken and instantons exist will be analyzed at length. The main 
motivation to discuss these ID SUSY QM models is to take profit of the knowledge acquired to study 
highly non-trivial 2D systems as those proposed in JH]- Another issue to be treated with care is the 
semiclassical behavior of supersymmetric quantum systems, this done with the help of the enlightening 
paper of A. Comtet et al. jllj . 
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2 Role of the Hamilton-Jacobi, Riccati and Poisson 
equations in SUSY quantum mechanics 



Let us start with a natural Lagrangian of one degree of freedom and the action functional: 



S 



dt 



mdx dx . 

r-—, V(x;X,k) 

2 dt dt K ' ' ' 



[A] = ML" 2 T~ 2 , [k] = MT~ 



(1) 



We shall consider potential energies V(x; A, k) that depend on two parameters A and k of dimensions 
given in and we shall introduce the non-dimensional variables: x — > y^-x, t — > y^r-i, V(x; A, k) = 
j^V(x), such that the action and the Hamiltonian read (non-dimensional variables will be used in 
what follows): 



S 



i 5 - 

k2m,2 



dt 



1 dx dx 
2~dt~dt 



V(x) 



H(p, x) 



A 



p 2 + V(x) 



P 



dL 

dx 



dx 
~dt 



2.1 One-dimensional J\f = 2 SUSY classical mechanics 



A M = 2 supersymmetric extension of a classical mechanical system of one degree of freedom is 
constructed as follows: 

1. We add two "fermionic" degrees of freedom to the "bosonic" degree of freedom with the real 
coordinate x. The fermionic coordinates form a Grassman Majorana spinor: 




>aPji 



+ o P e a = o 



Va,j8 = 1,2 



2. A superPoisson structure is defined in the phase superspace with coordinates p,x,0i,02 ■ Given 
two superfunctions F and G on the superspace, the Poisson superbracket 

{F G}p = _ + i y F JLJL G 

' dp dx dx dp ^—f d6 a d0 a 

is read from the basic brackets: Va, f3 = 1,2,, {p, x}p = 1, {x, x]p = 0, {p,p}p = 0, {0 a , 9p}p = i5 a p. 
Note that in the "soul" of the system - the subspace of the superspace spanned by the Grassman 
variables- the configuration space and the phase space coincide. The reason is that the Lagrangian 
ruling the dynamics of the fermionic variables is of first order in time derivatives. Thus, the time 
derivatives of Grassman variables will not appear in the Hamiltonian. 

3. The classical SUSY charges: Q\ = p9\ — ^&02, Q2 = P@2 + ^r6i, close the classical supersymmetric 
algebra: 

{Qi,Qi}p = {Q2,Q2}p = 2iH s 

4. The classical Hamiltonian H$ 

H, 



{Q a ,H s } P = 



{Qi,Q 2 }p = o 



IdWdW d 2 W 



1 

2^ 2 dx dx 



dx 2 



(2) 



is invariant by construction with respect to the super-transformations generated by Q\ and Q2- Besides 
the kinetic energy of the bosonic variables, there are two interaction energy terms in the Hamiltonian 
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([2"]) proportional to the (square of) the derivative and the second derivative of the arbitrary function 
W(x), usually referred to as the superpotential. 

Therefore, a given classical Hamiltonian: H = ^p 2 + V(x), admits an extension to a N = 2 
supersymmetric partner H$ if and only if the superpotential satisfies 

ldW dW 

Note that enters in H$ as the expectation value in Grassman states and disappears in a purely 
bosonic setting. 

Let us now consider the Hamiltonian for the "flipped" potential V{x) = —U{x) and the associated 
Hamilton-Jacobi equation: 

H F = -p 2 + U(x) ; -^+H F (—,x)=0 . 

The time-independence of the Hamiltonian suggests solutions of the form S(x,t) = —Et + W(x), 
leading to the reduced HJ equation: 

ldW dW rr . . 

— ; r- + U(x) = E 4 

2 dx dx K ' K ' 

Therefore, the superpotential is no more than the Hamilton characteristic function for E = of 
the mechanical system with flipped potential. In sum, to find the superpotential, allowing for the 
supersymmetric extension of a classical mechanical system, one must solve a related Hamilton-Jacobi 
equation, see Reference |7j. 

In general, for any E, the Hamilton characteristic function is: 

W(x;E)=± J dxy/2(E - U{x)) . (5) 

The energy E trajectories satisfy the ODE 

dx dW i—— — — , , f dx 

— = ±— = ±y/2(E - U(x)) => ± / — =t + t (6) 

dt dx V V K " J ^/2{E-U{x)) V ^ 

2.2 One-dimensional J\f = 2 SUSY quantum mechanics 

Canonical quantization of the above system to obtain the analogous M = 2 quantum supersymmetric 
system proceeds as follows, see, e.g., References [S], 0, HU and [T3] : 

1. Replace Poisson brackets by commutators for the bosonic variables and anticommutators for the 
fermionic variables: 

[x,p] = xp - px = ih , {0 a ,6p} = 9 a 0/3 + 6 a = -h5 al 3 

where the non-dimensional Planck constant K = — t^t has been introduced. 

2. We choose the coordinate representation for the bosonic variables but the classical Grassman 
variables become Fermi operators in the quantum domain: p = j-^, x = x, 0\ = Vhipi, 02 = vKip2, 
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The Fermi operators are represented on the Euclidean spinors in R 2 by the anti-Hermitian 2x2 Pauli 
matrices: 

V>1 = ^jfl , V>2 = , % = § Z = -\\ \ l J > {^1. ^} = > 1 § M = ihV3 ■ 

3. The quantum supercharges, Qi = -d\h£ - ih^, Q 2 = 2 h-^ + A^-, are 



and satisfy the quantum algebra: {Qi, Qi} = {Q2, Q2} = ^H s , {Qi, Q 2 } = 0, [Q\,H S ] = [Q 2 , H$] 
0, with the quantum SUSY Hamiltonian: 



\ _ 1 

U 1 ) ) ~ 2 



It is also interesting to work with non-hermitian supercharges Q± = \{Qi db iQ2), 

and reshuffle the quantum superalgebra in the form: {Q+, Q-} = 2hHs, [Q+, Hs] = [Q--> Hs] = 0. 

4. The quantum Hamiltonian is a block-diagonal 2x2 matrix differential operator h^ =0 ^ and h^ =1 ' 
are ordinary Schrodinger operators acting respectively on the subspaces of the Hilbert superspace 
labeled by the eigenvalues of the Fermi number operator: 

/=«-«+=*■ (° °) , h=^±>s 2 ) 

5. Wave functions in the subspaces with zero and one Fermi number annihilated respectively by 
Q + and Q_: Q+^\x) = 0, Q-^H X ) = °> 

are eieenfunctions of the Hamiltonian of zero energy. 

Therefore, 

= c( "*V W W ) , = C ( r A , „ 

W \ / ' W V expI-ft-^Ca;)] 

are the ground states of the supersymmetric quantum system if they are normalizable: dx e 2W ( x ^ 1 < 
+00 or 

J" R dxe~ 2W ( x )h 1 < +00. Note that either or 'I'q 1 ^ can be normalizable. 
2.3 The two-fold way to supersymmetric quantum mechanics 

Given a physical system, the issue of building the associated supersymmetric quantum mechanics can 
be addressed in two different ways. 
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• Quantization of a classical supersymmetric system. In the first method, it is assumed that the clas- 
sical supersymmetric extension has been performed. The identification of the classical superpotential 
requires that we must solve the ODE 

ldW dW 



V(x) 



2 dx dx 

the time-independent Hamilton-Jacobi equation © for zero energy and flipped potential energy. This 
idea has been applied to integrable but not separable systems with two degrees of freedom in Reference 
|19j . Canonical quantization, as in the previous Section, provides all the interactions in the quantum 
system 

V(0) M = IdWdW hfW n) IdWdW _ h£W 

1 1 2 dx dx ^ 2 dx 2 ' 1 ' 2 dx dx 2 dx 2 

in terms of the Hamilton characteristic function. 

• Supersymmetrization of a quantum system. The identification of the "quantum" superpotential 
would require one to solve one of the two Riccati differential equations 

IdWdW , hd 2 W Tr , , 

2^^ ±2^ = ^ ' (?) 
the sign marking the subspace where the the potential energy V is expected to act. There is no 
dependence on the Planck constant in the potential energy of any physically significant mechanical 
system. Therefore, we change the strategy and look for superpotentials that solve the Poisson equation: 

±^f(*) = V(x) , (8) 

dx* 

with the same criterion for the signs. Physically, this means that the Yukawa interactions provide the 
potential energy at stake. Mathematically, the solution of the Poisson equation (JHJ Wp provides a 
solution to a pair of related Riccati equations ©: 

i dw P dw P h<Pw P y(0) idWpdWp_hd 2 w P _ = ^ 1) 

2 dx dx 2 dx 2 2 dx dx 2 dx 2 { ' ' { ' 

for other related potential energies: V(o)(x), V^(x). Once again, the datum is V(x) in (JSJ) from which 
V^(x), V<V(x) are derived. 



3 Examples: Anharmonic oscillators of sixth-order 

To put these ideas to work, we choose as examples one-dimensional oscillators with terms proportional 
to x 4 and x 6 in the potential energy. Papers, reviews and even books dealing with the x 4 case abound. 
We shall discuss the x 6 case because it provides a splendid arena to disentangle two effects, instantons 
and spontaneous supersymmetry breaking, which in the x 4 case come together. The potential energies 
are: 

V{x-A,k) = ^x 2 (x 2 ±^j , V(x)= l -x\x 2 ±l) 2 , (10) 

describing respectively a single (+ sign) or triple (- sign) well. We shall only describe the first line 
of attack here from the solution to the HJ equation (where the potential energy is not found in the 
Yukawa interactions) and leave the Poisson route for another publication. 
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3.1 Quantization of classical supersymmetric sixth-order wells 
3.1.1 Single well 

1. Supersymmetric classical mechanics. The solution to the HJ equation for E = and U(x) = 
-\x 2 (x 2 + l) 2 is: 

W{x) = ±J dxx( X 2 + 1) = ± + yj 

The supersymmetric classical Hamiltonian and the supercharges read: 

H s = ^p 2 + -x 2 (x 2 + l) 2 Ti({Zx 2 + l)9 2 e l , Q a = p9 a =f x(x 2 + l)e a p6p 

In the "soul" of the related supersymmetric system with flipped potential, the Hamilton characteristic 
function and the trajectories are given analytically by hyperelliptic integrals: 

W(x;E)= f dx a/x 6 + 2x 4 + x 2 + 2E , f ^ =t + t 

J J Vx 6 + 2x 4 + x 2 + 2E 

For E = 0, there is only one constant trajectory, where the particle sits on the top of the potential: 

x(t) = 0, which is also the unique BPS trajectory of the supersymmetric classical system. 

2. Supersymmetric quantum mechanics. The quantum supercharges are: 

o n£ T x(x 2 + i)\ 6 = Jn( o -n£±x{* > + i) 

2 V h&±x{*> + l) ) ' ^ V 2 ^ - h £± X (^ + l) 



and the potential energies arising in Hs read: 

V®(x) = ^(x 2 (x 2 + l) 2 ±h(3x 2 + 1)) , vW(x) = ^{x 2 {x 2 + 1) 2 tH3x 2 + 1)) . (12) 
Thus, the zero energy ground states are: 

•f>(«)-c| ^{i^}) , <> w = cf ^ 

j \ exp{ T( t + '' } 

The supersymmetric quantum system has always one ground (BPS) state and supersymmetry is 
unbroken: if we choose W = ^ + ^ as the superpotential, the ground state belongs to the Fermi 
subspace - $g is not normalizable-, the choice of W = + forces a Bosonic ground state whereas 
^q 1 ^ becomes non-normalizable. 

One can guess the energy and the type of eigen-function of the next energetic states by looking at 
the "effective" potentials: 

V+(x) = ^(x 2 (x 2 + l) 2 -h(3x 2 + l)) , V^(x) = ^(x 2 (x 2 + l) 2 + h(3x 2 + l)) 



either or depending on the choice of W. The critical points of V±(x) are: xo = 0, xf 
d 2 V±, , oF d 2 V±, _ , - / , d 2 l/ : 



± 



2 W = lT3E , -^(xf)=4(l±3H+VT±m) , — ^ (xf ) = 4(1 ± 3ft - y/l ± 9ft) 
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V+ [x] 



V+[x] 



V+ [x] 



V+ [x] 






Figure 1: Potential energy V + and BPS ground state ^q(x) (red) for: 
(a) h = 0.001, (b) h = 0.1, (c) h = 2, (d) % = 4. 



xo is a minimum of V+ if K < ^ and becomes a maximum otherwise. x x are always imaginary roots 
but 

X<2 £1X6 real and become minima of V + for h > k. 



There is a unique minimum for V-, xq, and the wave function of the first level over the ground 
state is well approximated by a Gaussian around it: 



#5 (x) 







2 V 



\HttJ \ exp{-^x 2 } J 
The supersymmetric partner state in the subspace of ^(f is obtained by acting on with Q + : 



(13) 



(x) = Q+*Z{x) 



ujy | (x 3 + (1 + w_)x)exp[-^x 21 

47T 



(14) 



V-[x] 



V+ [x] 




V-[x] 




V+[x] 




V-[x] 




V+[x] 




Figure 2: Potential energies V±. Degenerate in energy , J , _ 1 (x) (green) 

E + _ _ _ 

and (x) (blue) wave functions: (a) h = 0.01, (b) h = 0.1, (c) h=l. 

4 2 
/ x j x \ _ 

3. Zero-energy ground state. The dependence of ^(f (x) = exp[ 4 2 ] on 7i is rather involved and 

can be described analytically through the asymptotic behavior when h — > h c and ^ c = is the classical 
value: 

2h 

exp|--(— +x')|< 1 + — • | — 



exp[-i( y + x 2 )] ~ Jim exp[-i(y + x 2 )] | 1 + i • 



+ 



1 - 2kr 



2h 



2h c V(^ + 



n c + 



i + m c (h c - 1) 

6M 



2h 



(f + * 2 ) 



M +••• 



It is also interesting to analyze how the norm of the BPS ground state ^^(x) depends on h: 



= / dxexp 



x 4 x 2 

2h ~ ~h 



dz exp 



hz 4 



X 



(15) 
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This non-gaussian integral is no more than the partition function 2(h) = N(h) of a QFT system in 
(0+0)-spacetime dimensions and Lagrangian |lUj : 

L = --if 2 - 4/ , z = ^= , X = 3h . (16) 

The partition function can be expressed as a series in h 



oo 



m = yj dz e~* , r dz z*» e~* = -^L - V5F , (17) 



m=0 



(4!) m m! ' (2m)\2 2m 

by performing infinite Gaussian integrals: 

2(h) _ « (-3^T . (4m)! = 1 _ 1 _ + 5^_ 2 _ 5^11 



(3^) + ^7(3^-^TF( 3/l ) +•- • ( 18 ) 



\/for ^ ( 4! ) m ^ ! (2m)!2 2m 8 V ; 3 • 2 7 v ; 3 • 2 10 

The expansion (|18|) of the partition function shows an essential singularity at h = -the classical 
limit- and it is an asymptotic series. The best approximation to the integral is reached by keeping a 
number of terms m,Q such that the quotient between two consecutive terms is of the order of one: 



(4mo + 3)(4mo + l) „1 
-\3h\ ~ hzrriQ w 1 =>• mo ~ — r 



4!(m + l) 11 u u 2ft 

and the error assumed by neglecting higher-order terms is bounded by exp[— ^r]. 

It is tempting to explain the pictorial description of the series using Feynman diagram technology. 
Writing the partition function in the form, 



; ^e(^ 2 )^ , (19) 



(-A)y? 4 v v (-AV 4 

v/A/3 ^7 m! 

one discovers the following Feynman rule: there is a single tetravalent vertex with a factor (—A). 
The lower-order terms in the series (|18|) correspond to the weights of the vacuum diagrams - the —A 
factor of the vertex divided by the combinatorial factor, the number of equivalent graphs of the same 
topological type- up to second order in perturbation theory shown in the next Table. 



Vacuum graph Weight Vacuum graph Weight 




8 



3.1.2 Triple well 



1. Super symmetric classical mechanics. The solution to the HJ equation for E = and U(x) 

2* 



]x 2 (x 2 - l) 2 is: 



W(x) 



dx x(x 2 



1) 



4 2 



4 2 

The superpotential is thus the "sombrero" potential. The supersymmetric classical Hamiltonian and 
the supercharges read: 



H S = \ P 2 + \x 2 {x 2 - I) 2 T i((3x 2 - 1)920! 



Q a = pOa T x(x 2 - l)e 



a/3 f/3 



Although feasible, we shall not attempt to search for trajectories with non-null Grassman degrees of 
freedom. 

It is interesting, instead, to look at solutions in the "body" of the related supersymmetric system 
with flipped potential because of their role in the quantum Hg system. The Hamilton characteristic 
function and the trajectories are given analytically by hyperelliptic integrals: 

dx 



W(x;E) = J dx \/x 6 - 2x 4 + x 2 + 2E , J 



t + to 



Vx 6 - 2x 4 + x 2 + 2E 

For E = 0, the integrations are easily performed and two kinds of trajectories are found: 

• Constant trajectories, where the particle sits on the top of the potential: x(t) = 0, x(t) = ±1. 

• Trajectories where the particle starts from a maximum of the potential at t = — oo and slowly moves 
to reach x = ±oo (infinite action) or another maximum (finite action) at t = oo. 

1 , . . . 1 



x 2 > 1 : x 2 (t) 



x < 1, instanton : x(t) 



1 _ e ±2(t+i ) ' ~ ~' • ~w ^ - e±2{t+to) 

The constant trajectories are special due to the fact that they are also zero energy (BPS) classical 



/VK7T 



Figure 3: (a) Potential energy U(x). (b) Hamilton characteristic function (superpotential) W ± (x). 
(c) Zero-energy, finite action, trajectories (instantons). 

solutions to Hs because the classical supercharges Q± = (p =F ix{x 2 — 1))0± are annihilated by them 
for any value of 6±. 

2. Supersymmetric quantum mechanics. The quantum supercharges are: 



2 V h-f±x(x 2 - 1) / V 2\Jif±x{x 2 -l) 
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and the potential energies arising in H$ read: 



V^(x) = ^(x 2 (x 2 -l) 2 ±h(3x 2 -l)) 
Thus, the zero-energy ground states are: 



VW( X ) = ^{x 2 {x 2 - l) 2 T K3x 2 - 1)) 



(21) 



^\x)=C 



exp{± 



n 







exp{ Mfci) } 



The supersymmetric quantum system always has one ground (BPS) state and supersymmetry is 

4 2 

unbroken: if we choose W = \ — ^- as superpotential the ground state belongs to the Fermi subspace 
,(0) 



is not normalizable-, the choice of W 



| — ^ forces a bosonic ground state whereas 



becomes non normalizable. 

Nevertheless, despite unbroken supersymmetry this system has instantons. To analyze the coexis- 
tence of these two phenomena one needs to study how 
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V + (x) = -(x 2 (x 2 -iy-h(3x'-l)) 



Vl(x) = -{x 2 {x 2 - l) 2 + h{3x 2 - 1)) 



2 V v " v ' 2 

evolve in response to changes in K. Note that either V+ or V- are either or V^ 1 ', depending on 
the choice of W. The critical points of V±(x) are: xq = 0, xf — i--*/ 2 v^l M^ ^ — / 2+y/i ±9h 



^±( Xo ) = l T 3H ; ^(xf) = 4(1 ±3/1- V^±9l) , ^±(xf)=4(l±3h+Vl±9ti) 



dx 2 



dx 2 



v+[x] 




V+[x] 



V+[x] 




V + [x] 




V+[x] 




Figure 4: Potential energy V+ and BPS wave function ^q(x) (red) plotted as functions of x 
for several values of H: (a) h = 0.001, (b) K = 0.1, (c) % = 2. 

x^ are always minima of V+(x), xq is a minimum of V+ if ^. < | but becomes a maximum if H > | 



3' 



and x^ are maxima of V+ for /i < i, not anymore critical point for h > i, see Figure 2. Therefore, 



because V+(a?o) 



(„±\ - 2+y / l+9l 
A x 2 J — 3 



2fr 



+ t>, is a false vacuum that decays to 



the true vacua x^ when h < |. The decay amplitude can be computed from the classical bounce for 
the flipped potential, starting and ending at xq, which is very well approximated by an instanton- 
anti-instanton configuration for small values of h. It is remarkable how well this behavior is described 
by the ground state wave function ^q(x); even more remarkable, ^q(x) also matches the expected 
behavior for h > g where there is no tunnel effect at all, see again Figure 2. 
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xq, however, is the absolute minimum of V-(x); if H < |, xf are also minima of V-(x), but 



V-(x ) 



4 < V-ixi 



2+a/1-9?1 
3 



-y/l-9h 



+ 2H 



i. If h > 5 xq is the single critical point 



(minimum) of V- (x). Therefore, the eigenfunction of the lowest eigenvalue of the Schrodinger operator 
with potential energy V-(x) is approximately a Gaussian centered at xq = 0: 



^" (x) "0 3 - exp{ -^ x2} ' 



(22) 



(x), the first eigenfunction of outside the kernel, lives in the subspace orthogonal to the sub- 
space of ty G (x). For h < grows from the decay of the false vacua x^ ruled by instantons/anti- 
instantons now starting and ending at x^. Mathematica drawings of these wave functions are offered 
in Figure 3. 



V-[x] 



lAAl 



V-[x] 



V-[x] 



V-[x] 





Figure 5: Potential energy V-(x) and wave function (x) (green) plotted as functions of x 
for several values of H: (a) h = 0.001, (b) h = 0.1, (c) h = 1. 

Acting on (x) with the supercharge operator Q + , an approximate eigenfunction of H$ is 



obtained in the subspace of ^ G (x). The supersymmetric partner of v I r ^ 1 (x) is thus 



*J (x) = Q+*f 1 (x) 



4-7T 



(x 3 — (1 — a;_)x)exp[- 



2h' 



-X 



E+ = E- X , (23) 



and Ef is the lowestdying eigenvalue in the subspace of the zero mode (ground state) . Plots of these 
"odd" wave functions are shown in Figure 4 for several values of h. The wave function has a node at 
the origin. 



on h is somewhat involved 



3. Zero- energy ground state. The dependence of (x) = exp[ 

and can be described analytically through the asymptotic behavior when h — > h c and h c = is the 
classical value: 



2h 



exp[ =( x 2 )l ~ lim expl — ■=-( x 2 )] < 1 + =^ • ( 

11 2h y 2 n n c ^o 11 h c y 2 n \ ^ hi \{^-x< 



+ 



l - 2h c 

2M 



2h 



+ 



i + m c (h c - 1) 

6M 



2h 



(f -* 2 ) 



+ ••• 



n 



V+[x] 



V+[x] 




E + 

Figure 6: Potential energy V+(x) and wave function (x) (blue) as a function of x 
for several values of h: (a) h = 0.001, (b) fi = 0.1, (c) K = 1. 



The norm of the BPS ground state ^q(x) is again a non-Gaussian integral. Denoting z = 2z = ip 



and |h = j3, we obtain: 



j dxexp + = e'iy/hj ebexpj--^ 2\^hz 3 - 2z 2 ^ 



(24) 



iV(^) = -2(fo) is the partition function for the Euclidean A(yj)o-model with spontaneous x — ► — x 
symmetry breaking in (0+0)-space time dimensions and the Lagrangian: 



1 







L = — ip" — —ip' 

v 4r 



3 ^3 



3! 



(25) 



Performing infinite Gaussian integrals 



EE 



e s ^4/3/3 m=0 2k=o 
one obtains the asymptotic ^-expansion: 



\2k- 



(4\) m m\ (3l) 2k (2k)l 



EE 



(-|9) m (-V5?) 



2A- 



(4m + 6A;)! 



efy/^rTs ^ o fr o (40 m m! (3!)2*(2fc)! (2m + 3fc)!2*»** 



(26) 



1 /? 5-7 o2 5-7-11 3 3- 5 -7 a2 5 a 5-7-11 2 5 • 7 • 11 • 13 • 17 » 
8 + 2^3 ~ 3~ + ^ 2 + -- + ^ + —2^f + 2^ ^ + 



Again, the optimum value of the number of terms of k type can be estimated. Keeping a fixed but finite 
value of m = mo such that m$ « ko, the quotient between two consecutive k = ko and k = ko + 1 
terms must be of the order of one: 



a mo+ko+l 



1 



(2& + 2)(2£;o + l) 
27f3k « 1 => k w 



(4m + 6fc + 6)(4m + 6fc + 5) • • • (4m + 6fc + 1) 
(2m + 3k + 3)(2m + 3k + 2)(2m + 3fc + 1) 



(3!) 2 2 3 



27/3 



and the error assumed by neglecting higher-order terms is bounded by exp[— 
Writing the partition function in the form 



2/3 

e 3 y/ 4/3/3 m=0 fc=o 



(~/3)y 4 x ... x (~/?)y 4 ' (-y^ 3 
y~v y~v / (4!) 



X • • • X 



3! 



771! 



(2*0! 



e-Vd^ (27) 
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one sees that the Feynman rules encompass one tetra-valent vertex and one trivalent vertex that are 



proportional respectively to (— (5) and — (\/3j3). Four-leg vertices come from 



(-i)y 4 

(4!) 



in the integrand 



of lj27l): three-leg vertices are due to 



-31 terms in (j27j) and only contribute in pairs. Comparison 

with the ^-expansion (|26[) shows that pictures of the k = terms, collected in the first two blocks of 
the first row, are provided by the diagrams shown in Table 1. Diagrams with one tetra-valent and 
two three-valent vertices, k = m = 1, shown in Table 2, provide the second block in the first row: 

7 -(3 2 . 



3/3 2 
"IB" 



3/3 2 

3~¥ 



W 2 
IF 



W 2 



W 2 



3/3 2 



3-5-7 , 
~2® 



Diagram 



Weight Diagram 



Weight 



3/^ 
" 2 6 



3/? 2 
3-2 5 



3^ 

' 2 4 

3^2 
' 2 3 



" 2 2 

3/3 2 
3-2 2 



In Table 3 only diagrams with tri-valent vertices, m = 0, are displayed: 



Vacuum graph 



Weight Vacuum graph 



Weight 



3p 



o 



3/3 
3-2^ 



CK) 



2o2 



¥3 



■o 



2o2 



3^/3 

JT2 5 



e 




3 2 /3 2 



2o2 



3 2 /3 
3-2 4 



QOO 



3 2 /3 2 
"2 s " 

3 2 /3 2 



3 2 /3 2 



3 2 /3 2 
"S - 2" 



Diagrams with two tri-valent vertices contribute: |f + = whereas the contribution of diagrams 



with four trivalent vertices is: 



+ 



3-2 5 



+ 



2 5 



+ 



~3~¥ 



+ 



3 2 /3 2 



+ 



2 5 



+ 



3^ 
3-2 



5-7-11 fl 2 



~2^ 



(3\ 
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4 Supersymmetric WKB approximation 



The semiclassical regime is characterized by the inequality: 



d 2 W 



dx 2 



(x) 



dx 



\n(x) ■ n(x)\' 



Thus, 



A 



« 



2V(x) 



dV/dx 



2n ^Wix) 

is satisfied in the limit of short wave lengths. To obtain the WKB eigen-functions of the SUSY 
Hamiltonian in, e.g. , the subspace for which the zero Fermi number is zero - because supersymmetry 
WKB eigenfunctions of non-zero energy in the Fermi sector are given automatically - one starts from 
the Wentzel-Krammers-Brillouin ansatz in the classically forbidden region E < V(x): 

W E (x) 



^(x,t) = A(x) ■ exp[- 
The Schrodinger equation for V+(cc) becomes 



z ' d 2 In A , , d In A , , d In A , , 
h 2 ( -r^-Or) + • -7—(x) 



e » 



(28) 



dx 2 



dx 



dx 



- H (^^(x) + 21nA(x) • ^r^(x) - ^(x)] + W E (x) ■ W E {x) - Sl(x) • + = (29) 
\ ax^ dx dx J 

with three terms of respectively order 2,1, and in h. The usual WKB strategy starts by solving the 
equation (|29j) for the ^.-independent terms to find: 



W E (x) = J dx y/n(x) ■ n(x) - 2E 



with the novelty with respect to the non SUSY case that the turning points are those corresponding to 
V(x), rather than those set by the effective potential V+(x). The second step is to plug this solution 
into the equation for the terms proportional to h: 



dlnA 

dx 



(x) 



1 



ft(x) 



y/tt{x)Q{x) - 2E n(x)n(x) - 2E 



dn 

dx 



(x) 



Integration of this equation provides the SUSY WKB wave functions: 



A{x) oc 



1 



|fi|(x) + y/n(x)£l(x) - 2E 



(30) 



(n(x)n(x) - 2E)4 

Note the other difference: in the non-SUSY case the numerator of this expression is 1. In the classical 
allowed regions, 2E > Q 2 (x), however, the WKB ansatz reads, 

\W{x)\. 



tp (x) = ^4(x)exp[— i- 



and one obtains: 



A{x) 



[2E-Q 2 (x)]i 



- ■ exp 



1 O(x) 
— arcsm 

2 J2E 



(31) 
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To match the WKB wave functions (|28|) and 1)3 1JI analytically at the classical turning points x = a < 
x = b, such that Q 2 (a) = 2E = 2 (6), the following supersymmetric quantization rule is required: 



6 i 

dx [2E - Q 2 (x)] 2 = nirh n £ Z+ . (32) 



The appearance of the numerator in ()3(J|) is magic: firstly, because this term modifies the process of 
analytic continuation necessary to match the exponential and periodic WKB wave functions at the 
turning points in such a way that the term that appears in the non SUSY version of 1)32)1 does not 
enter the SUSY case. To obtain the WKB wave function in the classically allowed region 



(2E-n 2 )± 1 J 



from the WKB wave functions in the forbidden regions 



Vn 2 -2E + \n\ lraj)/ -—-— „ Jvn 2 -2E + \n 



^e {x) = C TL ZL2. e -U: & VWPS . ^Er x) = ell ZJ_i. e -jU" ^'Vw 

(0 2 -2£)3 (^2-2^)4 

one chooses paths in the a>complex plane that goes around the turning points a and b at great distance, 
either in the upper or the lower half-planes. Unlike to the non-SUSY case, there is no e~ l ~± factor left 
and two wave functions are obtained in the classically allowed region, one from the left and the other 
from the right: 



\ J2E - n 2 (x) + in(x) i rx , 

^ E {x) = C^— j cos[i / dx' J2E - n 2 (x')] 

(2E-n 2 (x))* hJb 



Jj2E-n 2 (x)+in(x) i r a 

if> E (x) = C'±— 1 cos[i / dx' J2E-n 2 (x')] 

{2E-n 2 {x))\ KJx 

These expressions are identical if and only if 1)32)1 holds. Secondly, E = is a solution of 1)32)1 for 
n = 0, whereas l|2"H|l becomes the exponential of the superpotential: the exact ground state is a SUSY 
WKB wave function ! 

4.1 WKB analysis of the single well 

We shall consider as examples non-harmonic oscillators of fourth order to avoid hyperelliptic integrals 
and deal with (slightly!) manageable expressions. In the case of a single well with potential energy 
V(x) = ^(x 2 + j) 2 we have, using non-dimensional variables: 

V(x) = \{x A + 2x 2 + 1) , n(x)=x 2 + l , W{x) = ^- + x 
The turning points are the real roots of the quartic equation: 



x A + 2x 2 -a = , a = 2E-\ , x± = ±v/-l + y/l + a . (33) 
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The super symmetric quantization rule is therefore: 

f x + , 

I(E;x—,xj r )= I \/ a — x 4 — 2x 2 dx = nirh 



Denoting A± = ±1 + yjl + a, the definite integral in (|37|) reads: 



I(E;x-,x+) = -VT y/a + U- 



-4, 



(34) 



(35) 



where IEC(/c 2 ) and E(fc 2 ) are respectively the complete elliptic integrals of first and second type. This 
result is shown in Figure 8. 



V+[x] 



V+[x] 



V-[x] 



V-[x] 







Figure 7: Graphics of V ± {x) for H = 0.1 -(1),(3)- and h = 1 -(2), (4)- 




Figure 8: Mathematica plot of I(E; X j X _|_ ^ clS cl function of a and intersection with nirh for low n and 

ft = 1 (left) and h = 0.1 (right) 

The first three (double, see Figure 7) eigenvalues for h = 1 and ft = 0.1 are: -Ei = 2.18674^-, 
E 2 = 4.23942 ^, £ 3 = 6.5444 ^, and E 1 = 0.64500 ^, E 2 = 0.78289 ^, £ 3 = 0.95403 *f , respectively. 

4.2 WKB analysis of the double well 

For a non-harmonic oscillator of fourth order and a double well things are even more difficult. The 
potential energy is V(x) = ^(x 2 — -|) 2 , such that in non-dimensional variables we have: 

V(x) = hx 4 - 2x 2 + 1) , n(x)=x 2 -l , W{x) = ^--x 
The turning points are the real solutions of the quartic equation: 

x 4 -2x 2 -a = , = -JlTVl + a , x +T = a/ 1 T Vl + a ■ (36) 
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For a > there are only two real roots and the supersymmetric quantization rule reads: 

fX + + 



I(E;x- + ,x ++ ) 



\J a — x 4 + 2x 2 dx = mrh 



(37) 



The computation of I(E; x |_, x ++ ) is qualitatively identical to the previous case and results are shown 

in Fig. 10(left). 

If — 1 < a < things are more difficult: there are four turning points, four real roots, and the 
quantization rule splits into two equations: 



I(E;x ,x.\ )= / \J a — x 4 + 2x 2 dx = nnH = I(E; x 

J X 

The definite integrals in (|38|) now read: 



f x++ / 

/ v a — x 4 + 2x 2 dx 

J X |_ 



(38) 



I(E; x , x_| ) = I(E; x |_, x + _ 



2a 



-E 



-A- -Aa 



arcsm . 



A, 1 A. 



-A- -A. 



arcsm . 



A+ ' A. 



(39) 



Note that incomplete elliptic integrals of the first, W(u,m), and second, M.(u,m), type also enter. In 
any case, it is possible to plot these functions of a and find the intersection points determining the 
spectrum. 



V+ [X] 



V+[x] 



V-[x] 



V-[x] 







Figure 9: Graphics of V ± (x) for h = 0.1 -(1),(3)- and h = 1 -(2), (4)- 



10 12 14 



-1 -0.8 -0.6 -0.4 -0.2 



Figure 10: Mathematica plots of I(E; x ^,x ++ ) for h = 1 and a > (left) and I(E; x + _) for ft = 0.1 and 

— 1 < a < (right) as function of a. The intersection points with rmh giving the eigenvalues are also shown. 

The first three eigenvalues for a > and h = 1 are: E x = 0.82272 *f , E 2 = 2.08330 ^, E 3 = 
5.63830^. 
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In the case of — 1 < a < eigenvalues only exist if H < 0.95. Application of rule (|38j) for the 
turning points on the left gives: E x = 0.19183 E 2 = 0.36384 E 3 = 0.49993 

Because of formula (|39|) the choice of pair of turning points is irrelevant; E\, E 2 , E$, etcetera, are 
eigenvalues of the Schrodinger equation for both V+(x) and V-. 

5 Outlook 

The next step is to study physical systems of two degrees of freedom. It is tempting to start by 
discussing problems of this type in Hamilton separable systems. Following the works |14j-|15j on 
supersymmetric quantum mechanics in more than one dimensions, the general structure of supersym- 
metric classical and quantum Liouville systems has been described in References |13j and J2]. An 
important example of this kind of systems is the supersymmetric classical and quantum hydrogen atom 
respectively analyzed by Heumann ^\ and Kirchberg et al JH]- It seems also plausible to address 
similar issues in non-separable but integrable systems as those proposed in [T9*] . 

6 Acknowledgements 

We thank M. Ioffe for enlightening conversations on supersymmetric quantum mechanics and a crit- 
ical reading of the manuscript. This work has been partially supported by the Spanish MEC grant 
BFM2003-00936 and the project VA013C05 of the Junta de Castilla y Leon. 

References 

[1] J. F. Carihena, Motion on Lie Groups and Applications in Classical and Quantum Mechanics, 
Anales de Fisica, Monografias 6 (2002)211-226, Real Sociedad Espaola de Fisica 



[2] 
[3] 
[4] 
[5] 

[6] 
[7] 
[8] 
[9] 
[10] 



J. F. Carihena and J. Ramos, Reviews in Math. Phys. A12 (2000) 1279-1304 
J. F. Carihena, D. J. Fernandez, and J. Ramos, Ann. Phys. 292 (2001) 42 

F. Cooper, A. Khare and U. Sukhatme, Phys: Rep. 25(1995) 268 

G. Junker, Supersymmetric methods in quantum and statistical physics, Springer Verlag, Berlin, 
1996 



A. Wipf, Non-perturbative methods in supersymmetric theories, hep-th/0504180 



H. Goldstein, Classical Mechanics, Addison -Wesley, Reading, Massachusetts, 1980 

R. Casalbuoni, Nuovo Cimento 33A (1976)115 

N.V.Borisov, M.I.Eides, M.V.Ioffe Theor.Math.Phys. 29 (1976)906 

R. Borcherds and A. Barnard, Lectures on Quantum Field Theory, math-phys/0204014 



18 



[11] A. Comtet, D. Bandrauk and D. Campbell, Phys. Lett. B150 (1985) 159 

[12] A. Alonso, M.A. Gonzalez, M. de la Torre and J. Mateos Guilarte, J. Phys. A 37 (2004) 10323. 

[13] A. Alonso, M.A. Gonzalez, J. Mateos Guilarte and M. de la Torre, Ann. Phys. 308 (2003) 664. 

[14] A. A. Andrianov, N.V. Borisov. M.I. Eides and M. V. Ioffe, Phys. Lett A109 (1985) 143. 

[15] A. A. Andrianov, N.V. Borisov. M.I. Eides and M.V. Ioffe, Theor. Math. Phys. 61 (1984) 963. 

[16] P.F. Byrd and M.D. Friedman, Handbook of Elliptic Integrals for Engineers and Scientist, 
Springer- Ver lag 1971. 

[17] R. Heumann, J. Phys. A 35 (2002) 7437. 

[18] A. Kirchberg, J.D. Lange, P.A.G. Pisani and A. Wipf, Ann. Phys. 303 (2003) 359. 

[19] M.V. Ioffe, J. Mateos Guilarte, and P.A. Valinevich, Two-dimensional super symmetry: from susy 
quantum mechanics to classically integrable models, Ann. Phys., to appear, hep-th/0603006. 



19 



